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To study the ground state of ANNNI chain under transverse field as a function of frustration 
parameter k and field strength F, we present here two different perturbative analyses. In one, we 
consider the (known) ground state at k, = 0.5 and F = as the unperturbed state and treat an 
increase of the field from to F coupled with an increase of k from 0.5 to 0.5 + rT as perturbation. 
The first order perturbation correction to eigenvalue can be calculated exactly and we could conclude 
that there are only two phase transition lines emanating from the point k — 0.5, V = 0. In the 
second perturbation scheme, we consider the number of domains of length 1 as the perturbation 
and obtain the zero-th order eigenfunction for the perturbed ground state. From the longitudinal 
spin-spin correlation, we conclude that floating phase exists for small values of transverse field over 
the entire region intermediate between the ferromagnetic phase and antiphase. 

PACS numbers: 64.60.Fr, 75.10.Jm, 05.50.+q 
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I. INTRODUCTION 

The transverse Axial Next-Nearest Neighbour Ising 
(ANNNI) model is one of the simplest Ising models that 
contains tunable frustration and tunable quantum fluc- 
tuation. In one dimension, it is defined (for spin=i) by 
the Hamiltonian [l[ 
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)-r£. 



(i) 



Here the first term describes a ferromagnetic nearest- 
neighbour interaction of strength J(> 0) between lon- 
gitudinal components of spin s z (= ±1), the second term 
describes an antiferromagnetic second neighbour interac- 
tion of strength k J also in the longitudinal direction, and 
the third term describes an external field in the trans- 
verse direction of strength T. The ratio K (> 0) is called 
the frustration parameter. This Hamiltonian describes 
a classical ANNNI chain subjected to a transverse field 
as well as a transverse (nearest-neighbour) Ising model 
with an additional frustrated second nearest-neighbour 
interaction. In this paper, we shall deal with the ground 
state (zero-temperature) phase diagram of the Hamilto- 
nian Tt. It is known [J Q that for n = (only nearest- 
neighbour interaction), there is a ferromagnetic to para- 
magnetic second order phase transition at T = J. On the 
other hand, for T = (classical ANNNI chain) 0, H the 
ground state is ferromagnetic for k < 0.5 and antiphase 

(+ H type) for k > 0.5, with a 'multiphase point' 

at k — 0.5. At the multiphase point [3| the ground- 
state has very high (~ g , where g is the golden ra- 
tio (V5 + l)/2) de generacy as any combination of an- 
tiphase and ferromagnetic patches will serve as ground 
state configuration. From approximate analytic and nu- 
merical approaches [l|, 0], early studies had proposed a 
phase diagram (Fig. [T|) that consists of ferromagnetic, 



paramagnetic and antiphase regions, alongwith a region 
of what is called floating phase. The n th neighbour spin- 
spin correlation function in the longitudinal direction 
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decays exponentially with distance in the ferromagnetic, 
antiphase and paramagnetic regions but decays alge- 
braically in the floating phase. The evidence for the 
presence of floating phase is provided by quantum Monte 
Carlo simulation [5J] and by exact numerical diagonalisa- 
tion of small systems [3, || , one of which Q extends even 
upto a length of 32. 

Recently, a question has arisen on the existence of the 
floating phase, as a recent study [7[ has claimed that 
the floating phase exists only within a strip of infinites- 
imally small width. The transverse ANNNI chain is re- 
lated to the two-dimensional classical ANNNI model by 
Suzuki- Trotter transformation. The phase diagram for 
this model is also similar to Fig. [T] (with T replaced by 
temperature). Here also the earlier studies 0] support 
the existence of floating phase in view of Monte Carlo 
simulations by Selke and others [8| and approximate an- 
alytic calculations by Villain and Bak M, while a recent 
measurement of dynamical exponent |10| and a density 
matrix renormalisation group analysis |FJ| claim that the 
floating phase exists, if at all, along a line only. We have 
also studied [l2j some static and dynamical properties of 
the 2D ANNNI model by Monte Carlo simulation. Our 
study also indicates that the floating phase exists, if at 
all, only along a line. The question is to whether the 
floating phase exists only along a line or extends over a 
region in the case of transverse ANNNI chain is the basic 
motivation of this paper. 

In this paper we perform two perturbation calculations 
with an aim to have an idea about the phase that exists 
near the multiphase point for non-zero transverse field. 
Our conclusion is that there exists floating phase over a 
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FIG. 1: 



Schematic phase diagram of the transverse ANNNI model 
according to early investigations [1,2]. 
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FIG. 2: 

Schematic phase diagram of the transverse ANNNI model 
according to the present work. 



region extending from ferromagnetic phase to antiphase 
(for small values of T) and the phase diagram looks like 
Fig. [2] This result is in contradiction with previous re- 
sults as none of the previous studies had predicted float- 
ing phase for k < 0.5. We must mention that all our 
results are true only at "small" values of T and the quan- 
titative details of this diagram is not reliable at T ~ 1. 
However, the topological structure of the diagram should 
be correct. This article is organised as follows. In Sec- 
tions II and III, the first and the second perturbation 
schemes will be presented (respectively) and we shall con- 
clude with some discussions in Sec. IV. 



II. THE FIRST PERTURBATION SCHEME 

A. Principle 

We shall now present a first order perturbation calcu- 
lation around k = 0.5 by rewriting the Hamiltonian Ti of 
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FIG. 3: 

Scheme of the first perturbation treatment. The transverse 

field is altered from to a (small) value F and the 
frustration parameter k is altered from 0.5 to 0.5 + rr. The 
parameter r may vary from — oo to oo. 
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FIG. 4: 

One possible type of phase diagram of the transverse 
ANNNI model. 



Eq. (1) as 

7~t = H c i + Hp 

where 

«d = -JB^i-^ +J ] (3) 

j 

and 

«P=E r [- s i+^+J (4) 

j 

with r = J(k — 0.5)/r. We always take T to be positive, 
so that negative values of r imply k < 0.5. 

We shall treat the TC P as a perturbation over Ti c i- Ob- 
viously Hci is the longitudinal (classical) part of H at 
K = 0.5 and Ti. p includes the transverse part and an 
increase of k (from 0.5) by an amount rr. One must 
note that for a perturbative treatment T has to be small 
but r need not be so. In the (k, F) phase space, the 
perturbation takes us from the multiphase point (0.5, 
0) to a point (0.5 + rr, T). As r varies from —oo to 
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oo, the perturbation takes us from the ferromagnetic to 
the antiphase region (Fig. [3]). The first order correction 
(E^ 1 ') to ground state energy can be calculated exactly. 
This quantity is studied as a function of r, because if 
there is a non-analyticity at r — r$, then there will be a 
phase transition line emanating from the (0.5, 0) point 
at an angle tan _1 (l/ro) with the increasing k direction. 
It is found that E { • ) as a function of r is non-analytic 
only at r — —0.5 and 1. Therefore, only two lines of 
phase transition meet at the multiphase point. This im- 
plies either that the floating-paramagnetic transition line 
touches tangentially the floating-antiphase line at this 
point (Fig. [3} or that the floating phase extends over 
the entire region from antiphase to ferromagnetic region 
(Fig. It will be shown that a study of susceptibil- 
ity and mass-gap rules out the possibility of Fig. 0] and 
points to the validity of Fig. [21 Thus, the final conclusion 
from the first perturbation scheme is that the boundary 
between the ferromagnetic (antiphase) and the floating 
phase emanates from the multiphase point at an angle 
tan -1 —2 (tan -1 1) with the k axis. The second pertur- 
bation scheme also supports this conclusion and addition- 
ally enables us to obtain an approximate phase diagram. 
We shall now transform our problem to a certain prob- 
lem for the nearest neighbour transverse Ising model and 
then through a study of the nearest neighbour model, 
obtain the expression for E^\ 

B. Transformation to Nearest-neighbour Model 

We start by noting that at re = 0.5 (the 'multiphase 
point' (H) the ground state of H c i, the unperturbed 
Hamiltonian, is a highly degenerate state and any spin 
configuration that has no spin-domain of length unity 

(i.e. no+H h+or 1 type arrangement) 

can be the ground state. The number of domain walls 
is immaterial and can be anything between and N/2, 
N being the total number of spins. (Of course, for pe- 
riodic boundary there can be only an even number of 
walls.) Let us denote the set of all such configurations as 
S. Also, let the population of this set be v which inci- 
dentally is of the order of g N as mentioned in Sec. I [|| . 
Now, the first-order correction to the eigenvalue are the 
eigenvalues of the v x v matrix P, whose elements are 

P a p ee ( a \n p \p) 

where \a) and \/3) are configurations within S . The 
matrix P can be easily given a block diagonal structure. 
Note that sJ|/3) £ S if and only if the j-th spin lies at 
the boundary of a domain, and the domain too has length 
larger than 2. Also, in such a case, operating on \(3) 
will translate the wall at the left (right) of the j-th. site by 
one lattice spacing to the right (left). This immediately 
leads us to the important conclusion that P a p ^ if and 
only if | a) and \/3) have equal number of domain walls. 
Thus, we can break up S into subsets S(W), where S(W) 
contains all possible spin distributions with W walls (W 



= 2, 4, ■ • ■ N/2). Now, the u x v matrix P gets block- 
diagonalised into matrices of size i>w x ; 

P a p{W) ee (a\H p \p) (5) 

where vw is the population of S(W) and \a), \0) <G S(W). 
We now observe that the longitudinal term in Ti. p only 
contributes a diagonal term rT(N — AW) to P a p(W), so 
that one can write 

P(W) =M{W)+rT{N- AW)1. (6) 

Here 1 is the vw x v w unit matrix, M a p(W) = {a\TL q \(3) 
and Ti q = — T sj is the transverse part of 7i p . Thus 
the non-trivial problem is to solve the eigenproblem of 
M(W). 

To proceed with the matrices M(W), let us construct 
from each member \a) of S(W) a configuration \a!) by 
removing one spin from each domain. The total number 
of spins in \a') will obviously be N — W = N' , say. Such 
a transformation was also used by Villain and Bak Q 
for the case of two-dimensional ANNNI model. The set 
S'(W) composed of the states \a') is then nothing but 
the set of all possible distributions of N' spins with W 
walls, with no restriction on the domains of length unity. 
It is easily seen that the matrix 

M' ap {W) ee (a'\H q \(3') 

will then be identical with M(W) since the element 

N' 

(«'iE s ^'> 

3=1 

is non-zero when and only when 

N 
3=1 

is non-zero. The eigenproblem of M'(W) becomes simple 
once we we observe that S'(W) is nothing but the set of 
degenerate eigenstates of the Hamiltonian 

N' 
3=1 

corresponding to the eigenvalue 

E w = -J(N' -2W). (7) 

Thus, if we perturb H.' by TL q , then the first-order per- 
turbation matrix will assume a block diagonal form made 
up of the matrices M'{W) for all possible values of W. 

C. Study of the Nearest-neighbour Model and 
Further Analysis 

To solve the perturbation problem for H.' Q + TL q , we 
note that this Hamiltonian is the same as Tt of Eq. (1) 
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with k = 0, namely, 

i=i 

This is the Hamiltonian for the standard transverse Ising 
model. The exact solution for this Hamiltonian is readily 
available 0, [l3|, [l4[ ■ The exact expression for the energy 
eigenstates are (for periodic chain in the thermodynamic 
limit [H), 

E = 2Tj2^k (8) 

k 

where may be 0, ±1 and fc runs over N' /2 equispaced 
values in the interval to ir. Also, stands for a/(A 2 + 
2Acosfc + 1), where A is the ratio J/Y. For T = 0, the 
energy E must be the same as % of Eq.(7), so that 

2 = -2W). (9) 

Clearly, the different values of the quantity 

V 0,1 / r=o fc=0 

correspond to the first order perturbation corrections to 
the different levels. They are therefore also the eigenval- 
ues of the M' matrix. Thus the eigenvalues of the matrix 
P(W) of Eq.(6) are 

E P =rT(N- 4W)+2r^£ fe cosfc (10) 

k=0 

Keeping N fixed we have to find, for which value of W 
and for which distribution of Ep is minimum subject 
to the constraint (9). For a given value of ^£fc ; this 
minimisation is achieved if -1 values of accumulate 
near k = and +1 values near k = tt. Let the desired 
distribution be 

( -1 for k = to 9 
ik= I for k = 9 to <f> (11) 
[ 1 for k = (f> to 7r 

Eq.(9) now gives 

AT/JV' = (4tt - 9 - <jj) /2?r (12) 
and one obtains, 

at 

E P = r(47T- 36>- 30) + 2(sin0 + sine/.)] 

47T — — (f> 

This quantity attains a minimum value only when 9 and 
(/> are equal and their common value (<fo, say) satisfies 
the condition 

27rr = sin^o + (27T — o )cos0o, (13) 



The minimum value of Ep is given by, 

=-jVT[3r-2cos0 o ]. (14) 

This is the exact expression for the first order perturba- 
tion correction to ground state energy. It is easily seen 
that for r < —0.5, that is, for T/J < (1 — 2k), one has 
4>o = tt and W = (ferromagnetic phase), while for 
r > 1, that is T/J < (re — 0.5), one has </>o = and 
VF = 7V/2 (antiphase). As r varies from —0.5 to 1, <fio 
gradually changes from tt to according to Eq.(13). The 
second derivative of E^ with respect to r blows up at 
r = —0.5 and 1, indicating two critical lines there. One 
can see from Eq. (13) that except for these two values, r 
is an analytic function of 4>o and hence, by Eq. (14), E^ 
is also an analytic function of r. We have also checked 
explicitly (Appendix A) that no higher derivative of E^ 
with respect to r blows up at any other value of r. 

As mentioned in Sec. IIA, we can now conclude that 
the phase diagram is either like Fig. [2] or like Fig. |4j 

D. Study of Longitudinal Susceptibility 

We shall now show that an analysis of longitudinal 
susceptibility points to the possibility of Fig. O rather 
than of Fig. g] 

Let us call the eigenstate of Ti.' + H q corresponding 
to 8 = <f> = <fio as This state will be composed of 

the spin-distributions that belong to S'{W) and can be 
written as 

where \f) runs over all the states in S'(W). If we con- 
struct from each state \ another state \j) by adding one 
spin to each domain, and then combine these states with 
the same coefficients, then we arrive at a state Y) j a,j \j) 
where a,j — ay . This is an eigenstate of M(W) and hence 
of P(W) (see Eq. (6)) and this eigenstate is nothing but 
the zero-th order eigenfunction IV' ) for the perturbed 
ground state of H c i + T~t P - One should observe that al- 
though the spin-spin correlation may not be equal for |j) 
and \j'), the longitudinal magnetisation M z must be the 
same for them (- equal number of positive and negative 
spins are added while transforming \j) to \j')). Thus, the 
longitudinal susceptibility 

Xz ex (Ml) ~ (M z ) 2 

of |-0' ^) must be the same as that of The spin-spin 
correlation 

C z {n) = {stsi +n ) 

for I?//} may be calculated (see Appendix B). In the case 
of r < J, for the entire range < (j>o < tt, the correlation 
is 

C z (n) = A-^= cos[n(7r - O )] 
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where A is a constant. This is clearly a floating phase. 
The susceptibility \z is hence infinity for both the states 
\tp') and This leads us to the conclusion that the 

zero-th order eigenstate is in floating phase and hence, at 
least for small values of T, the ground state of transverse 
ANNNI chain must be a floating phase for all values of 
r between —0.5 and 1. Of course, for large values of T 
the perturbation corrections may cancel the divergence 
of susceptibility and lead to a paramagnetic state. 

One must note that it is difficult to derive an exact re- 
lationship between the correlation in state |^°)) and the 
same in state | tp'}. Although a similar relationship was 
obtained by Villain and Bak we do not extend that 
derivation here. It is interesting to note that Villain and 
Bak (and also Uimin and Rieger [J]) assumed the wave 
number to be equal to the number of domain walls per 
site, but we shall soon find that this conclusion disagrees 
with the results obtained from our second perturbation 
scheme. 



E. Study of Mass Gap 

We shall now show that an analysis of mass gap also 
points to the possibility of Fig. [2l rather than of Fig. [4] 
thus agreeing with the conclusion from the analysis of 
longitudinal susceptibility. 

One signature of floating phase or diverging correlation 
length is vanishing mass-gap 0, [l6| . Let us now study 
the first order (in T) correction to the mass-gap. Since 
the first order correction to all energy states is given by 
the different possible values of Ep of Eq. (10), the correc- 
tion to the energy of the first excited state is the smallest 
possible value of Ep apart from the ground state E^\ 
To find the lowest excitation over the ground state, we 
note that such excitation is possible either (i) by keeping 

£fc fixed and rearranging the values; or (ii) by al- 
tering 9 and (f> and thus altering ^ For (i) the lowest 
excitation will correspond to an interchange of +1 and -1 
at k = 4>o, which will lead to a mass gap (for the whole 
system, not per site) 



8?rrA sin < 



N'A 



4>o 



For (ii) the mass gap is : 



A (D 



d 2 E F 

ae 2 



Here 69 is the smallest possible deviation in 9 at <fro- As 
the smallest possible change in W, and hence in N' is 2, 
Eq.(12) (with 9 = <f>) tells us that 



2(2tt- 



ttN 



1 

N' 



This shows that for both the mechanisms (i) and (ii), 
the mass gap vanishes as N — > oo for all values of 



<fio between and n. This shows that for all values of 
r between —0.5 and 1 there must be floating phase for 
small T. 

For transverse Ising model (k = 0) the phase transition 
occurs at T — J and at T = J + e the first order (in e) 
correction to the mass gap is just T, and is thus nonva- 
nishing, indicating that the divergent correlation length 
does not extend beyond r — J. 



III. THE SECOND PERTURBATION SCHEME 
A. Principle 

In the previous section we have found that the phase 
diagram looks like Fig. [51 In order to reconfirm this re- 
sult explicitly, we have to know the eigenfunction with 
first- or even zero-th order correction. As it is difficult to 
calculate the eigenfunction for the perturbation scheme 
discussed above in Sec. II, we shall perform the second 
set of perturbation calculation now. 

Let us write the Hamiltonian Ti of Eq. (1) for k < 0.5 

as 



U = N(h + hi - Jk) 



with 



h 



-(1 



3 



N ^ 3 

3 



and 



hi 



1 



3 



s z jS z j+1 )(l-s 



(15) 



(16) 



(17) 



Clearly, ho represents the Hamiltonian of standard trans- 
verse Ising model with nearest-neighbour interaction and 
hi represents the number of domains of length 1. We 
shall treat the operator hi as perturbation on the Hamil- 
tonian ho- All the eigenstates of ho are precisely known 
[H[ and for each eigenstate one can readily calculate 
the expectation value of ho + hi and identify the eigen- 
state for which < ho + hi > is the lowest. This state 
is the eigenstate with zero-th order perturbation correc- 
tion. Therefore, the basic idea is to identify the ground 
state upto first order perturbation correction to eigen- 
value. We carry on this scheme for all values of k (< 0.5) 
and r and in each case calculate the correlation func- 
tion (as defined in Eq. (2)) and characterise the phase 
therefrom. 

Indeed, our results are reliable only for small values 
of < hi >. Hence, this quantity must be small for the 
low-lying eigenstates of h . As mentioned already, the 
operator hi basically counts the number of domains of 
length 1. For k < 0.5 the ground state of ho at T = is 
ferromagnetic and the low-energy states may be expected 
to have small < hi > and the perturbation treatment will 
be justified for small values of T. However, for n > 0.5, 
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the ground state of ho will be antiferromagnetic (having 
domains of length 1 only) and such perturbation scheme 
is not valid. 

For k > 0.5 one can break up TL as 



H = N(h' + h[ - 0.5 J) 
with the unperturbed Hamiltonian defined as 



4^ 



and the perturbation as 
1 



h[ 



j 



1 



-i)(l 



(18) 



(19) 



(20) 



The Hamiltonian h' represents a chain with only next- 
nearest neighbour interaction and the operator h[ counts 
(as for k < 0.5) the number of domains of length 1. The 
ground state of h' Q at T — is antiphase, and hence has 
< h[ >= 0, so that the perturbation scheme is justified. 
Following the principle for the previous case, one can 
identify numerically the eigenstate (of h' ) for which < 



h' 



h[ > is a minimum. 



From the correlation function 



of this state, one can characterise the phase. 

The phase diagram obtained from the analysis of the 
present section is shown in Fig. [5j It has the following 
features : (i) The floating phase exists and extends from 
the ferromagnetic phase to the antiphase for small T; thus 
Fig. [2] rather than Fig.[4]is the type of the phase diagram, 
(ii) There is a line along which transition from floating 
phase to paramagnetic phase takes place; this line for 
k < 0.5 (obtained using Eq.(15)) continues smoothly to 
the same for k > 0.5 (obtained using Eq.(18)). 

We shall now present the details of the calculation for 
the two ranges of k. 



B. For k < 0.5 

The Hamiltonian ho of Eq.(16) corresponds to nearest- 
neighbour transverse Ising model with (ferromagnetic) 
interaction strength (1 — 2k)J. Following Eq.(8), the 
eigenstates of this operator are 

1 



< ha >= r 



N/2 



k 



& Ah 



(21) 



with A fc = ^/(A 2 + 2Acosfc + 1) as in Sec. II, but A = 
(1 — 2k) J/T. The line corresponding to a constant value 
of A = is now a straight line at an angle tan -1 (2/ A) 
with the k axis (Fig. [|5]). For each eigenstate, the n-th 
neighbour longitudinal correlation function (defined by 
Eq. (2)) can be expressed as the Toeplitz determinant 

M 



C z {n) 



Gi 
G 2 



G_i G_2 
Go G_i 
Gi Go 



G n -1 G„-2 Gn-2 



G-n 
G-n 
G-n 

Go 



(22) 



where the elements are given by 

G i = -Jr E f^cos^ - k) + Acos(fcj)] 

k=0 k 



(23) 



Here, the wave vector k runs over the N/2 equispaced 
values in the interval to it and 6 can be or ±1 for 
each of the N/2 values of k between and ir.) 

Now, for any eigenstate of ho, the expectation value of 
hi can be calculated first by rewriting it as, 



1 = «jl£(l-2, 



and then evaluating the first and second neighbour cor- 
relation functions using Eq.(22). The result is, 



< hi >= kJ[(1 - G ) 2 - GiG_i]. 



(24) 



To find which distribution of gives the smallest value 
of < ho + hi >, we note that Go will be closest to 1, when 
the (—1) values of 6 accumulate near k = and (+1) 
values near k = ir. Although Eq.(24) does not clearly 
indicate that such a distribution of 6 will also lead to 
the largest values of Gi and G_i, we have verified by 
going through all the 2 N states (for N = 12) that indeed 
such a distribution leads to the lowest value of < ho + 
hi >. We can now assume the distribution of values 
to be the same as described in Eq. (11) and calculate 
< ho + hi > for given values of 9 and <p. For every 
k and r, we take a system of 1000 spins and consider 
all possible choices of 9, <fi and note the values (say, 9o 
and 4>o) for which < ho + hi > attains a minimum. For 
the entire parameter range investigated in this work, we 
found #o = 4>o- Eq.(ll) now reduces to (with = <j) = 4>o) 
a distribution, namely 



6 = 



-1 
1 



for k 
for k - 



to 4> 

4>o tO 7T. 



(25) 



This corresponds to our approximate ground-state eigen- 
function for 7i. We have also found that for this eigen- 
function, < hi > is non-zero. As < hi > is zero for all the 
ground state eigenfunctions at T = 0, we may conclude 
that the state we have determined as the ground state of 
ho + hi, is an excited state of ho- In this connection, we 
point out that for 2D ANNNI model, it has been proved 
that 17f if one neglects completely the domains of length 
1 (i.e. assumes < hi > to be zero), then one misses the 
floating phase, indicating that the floating phase consists 
of states that do not belong to the ground state at the 
(T, k) = (0, 0.5) point. 

For every value of k and T we now have determined 
the value of 4>o and Eq.(25) therefore gives us the ground 
state eigenfunction. This leads to an expression for the 
correlation function C z (n) of Eq.(22) in the form of a 
determinant, an analytic expression of which has been 
given in Appendix B for large values of n. The final result 
is as follows. For A > 1, that is, T/J < (1 — 2k), the value 
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of </>o turns out to be tt, leading to a ferromagnetic phase. 
For A < 1, that is, T/J> (1 — 2k), there are two regions. 
In one region </>o < tt and according to Appendix B, the 
correlation is 

C z {n)=A{^)±= 

indicating a floating phase with index 0.5 and no modula- 
tion. The parameter (f>o varies continuously as a function 
of k and r (Fig. [7]), but this variation affects only the 
amplitude A. In the other region, c/)q — tt, and the cor- 
relation is exponentially decaying, indicating a paramag- 
netic phase. The values of A and 4>q and the correspond- 
ing phases in different portions of the phase diagram is 
indicated in Fig. [H The phase diagram thus obtained 
agrees with the conclusions from our first perturbation 
scheme and is presented in Fig. [5] 

C. For k > 0.5 

The Hamiltonian h' of Eq.(18) corresponds to a chain 
(C, say) with only next-nearest neighbour interaction 
and can be broken into two independent transverse Ising 
chains (Ci and C2, say) each having nearest-neighbour 
antiferromagnetic interaction. Hence, the correlation 
J2 s j s j+i m C wm be zero and J^ . SjSj +2 in C will be 
the same as the nearest-neighbour correlation in the con- 
stituent chains C\ and C2. (Indeed, C\ and C2 will be 
uncorrelated because, any configuration of C\ will cou- 
ple with equal probability to two configurations of C2, 
one of which can be obtained from the other by reversing 
all spins.) Moreover, since each of the antiferromagnetic 
chains can be transformed to a ferromagnetic chain (Co, 
say) by simply reversing the alternate spins, the eigen- 
states of h' in C are related to the eigenstates of a nearest- 
neighbour ferromagnetic transverse Ising chain Co- We 
have then, 

<^>=r-^£> A * (26) 

and 

< h[ >= 0.5,7(1 - G ) (27) 

with Go defined by Eq. (23) but now with A = (k — 
0.5)J/r. As before, k runs over N/2 equispaced values 
in the interval to tt. The A = constant line is now 
a straight line at an angle tan _1 (l/A) with the k axis 
(Fig. O. As in the previous subsection, we look for the 
eigenstate, (i.e. the distribution) for which < h'o+h'i > 
will be lowest and observe that Go will be largest when 
the (—1) values of accumulate near k = and (+1) 
values near k = tt. Also, investigating all the states for a 
chain of only 12 spins, we find that the minimum value 
of < h' + h[ > corresponds to the distribution with no 
zero values of Thus, assuming the distribution of Eq. 



(25), we can find numerically the value of 9 (say, <po) for 
which < h' + h{ > will be lowest for a chain of 1000 
spins. It is found that this minimum correponds to a 
non-zero value of < h[ >, indicating that as for k < 0.5, 
the unperturbed state corresponding to the ground state 
is an excited state of h' . 

As in the previous subsection, we can use the expres- 
sion for C z (n) mentioned in the Appendix B as the cor- 
relation in the chain Co- The corresponding spin-spin 
correlation in the ANNNI chain C will be 

„ z , _ J for odd n 

WlW - I (_i)»C z (n/2) for even n 

In this expression, the (—1)™ factor takes account of the 
mapping between the ferromagnetic and antiferromag- 
netic chains. The final result is as follows. For A > 1, 
that is, r/J < (k — 0.5), the value of 4>o turns out to be 
tt, so that the correlation in the chain Co is constant and 
the correlation C A (n) corresponds to a perfect antiphase. 
For A < 1, that is, T/J> (n — 0.5), there are two regions. 
In one region <fio < tt and according to Appendix B, the 
correlation is 

C* A (n) = A(M cos (28) 

indicating a floating phase with index 0.5 and modulation 
q = tt/2. 

The parameter 0o varies continuously as a function of k 
and r (Fig. [5]), but this variation affects only the ampli- 
tude A. In the other region, <pQ — tt, and the correlation 
in Co is exponentially decaying, indicating a paramag- 
netic phase with the same modulation. The values of A 
and (fio and the corresponding phases in different portions 
of the phase diagram is indicated in Fig. [8] and the re- 
sulting phase diagram is presented in Fig. which agrees 
with the conclusions from our first perturbation scheme. 

IV. DISCUSSION 

In this section we shall point out some features of the 
phase diagram we have obtained. 

(1) In the floating phase (A < 1) the spin-spin correla- 
tion decays as 

C z ( n ) = A— 

with 77 = 0.5. It is interesting to note that on the bound- 
ary between the floating and the ordered (ferromagnetic 
or antiphase) state A is 1 and <p is or tt and the cor- 
relation decays algebraically with 77 = 0.25 as shown by 
Pfeuty Thus, the value of the index undergoes a 

non-analytic change at the boundary. 

(2) In the floating phase, the correlation is non- 
oscillatory (q = 0) for k < 0.5 and oscillatory (q = tt/2) 



8 



0.8 



0.6 



0.4 



0.2 



- Vk Paramagnetic > 

- V V / 


/ 

/ 


\ w 

\ Floating / 


Antiphase 


- Ferro V + 




, , v. , , 




3 0.2 0.4 0.6 0.8 1 


1.2 1.4 1.6 


K 




FIG. 5: 





17 J 



Phase diagram of the transverse ANNNI model as obtained 
from the first and the second perturbation scheme. The 
thick lines are exact results from the first perturbation 
scheme. The thin lines are approximate results from the 
second perturbation scheme. The thin line is not drawn in 
the region where the latter scheme becomes unreliable. 



X = const 




X = const 



FIG. 6: 

Plot of A = constant line with the k axis. 
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FIG. 7: 

Plot of <f>o vs. F for k < 0.5. Here A = (1 — 2k)J/T. 
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FIG. 8: 

Schematic phase diagram showing the values of A and (f>o in 
different portions. 
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FIG. 9: 

Plot of 0o vs. r for k > 0.5. Here A = (« - 0.5) J/F 



for k > 0.5. Hence, the line k = 0.5 is the "disorder 
line" [H across which the modulation wavevector suffers 
a sudden change from to tt/2. We have also mea- 
sured the perturbed energy < ho > + < hi > and 
< h' Q > + < h[ > on two sides of the disorder line 
for a given T (Fig. 1 10[) . This energy remains continuous 
but suffers a change of slope at k = 0.5. Since the first 
perturbation scheme does not indicate any phase transi- 
tion around n = 0.5, we guess that this change of slope is 
not the signature of any serious non-analytic behaviour 
but is only an outcome of the approximation inherent in 
perturbation calculation. 

(3) Within the floating phase, for k > 0.5, the nearest- 
neighbour correlation is zero (Eq. (28)) and the wave- 
vector is q = tt/2, implying that the wave number q/27T 
is equal to the number of domain walls per site. On the 
other hand, for k < 0.5, the nearest-neighbour correla- 
tion is not at all zero (Fig. ITT]) . but the correlation is 
non-oscillatory (q = 0), implying that the wave num- 
ber q/2n is not equal to the number of domain walls per 
site. One must note that (as mentioned above), the an- 
alytic treatments of Villain and Bak [9f, and of Uimin 
and Rieger Q assumes the equality of the wave-number 
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FIG. 10: 

Plot of ground state energy vs. k for P = 0.8 
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FIG. 11: 

Plot of nearest neighbour correlation as a function of V for 
A=0.75, 0.90 and 0.95 



and the number of domain walls per site. This assump- 
tion thus agrees with our results for k > 0.5 but not for 
K < 0.5. 

(4) The second perturbative scheme would be reliable 
only when the ratio < hi > / < ho > or < h[ > / < h' > 
(accordingly as n is < 0.5 or > 0.5) is small. We have 
included in Fig. O only those points where this ratio is 
less than some arbitrarily chosen number 1/3. This cri- 
terion made us unable to locate the boundary between 
the paramagnetic and floating phase near k = 0.5, T = 
resulting in a gap there in Fig. [5j Some better approxi- 
mation is needed to bridge this gap. 

(5) The first and second perturbation schemes can be 
compared and reconciled in the following way. For the 
first perturbation scheme, the unperturbed Hamiltonian 
corresponds to the k = 0.5, r = point, and the exact 
expressions for the first order correction to energy gives 
exactly the directions along which the boundaries of fer- 
romagnetic phase and antiphase emanate from this point 
(see Fig. 5). On the other hand, for the second perturba- 
tion scheme, the unperturbed Hamiltonian corresponds 
to the boundary lines of the ferromagnetic phase and the 
antiphase since on these lines < hi > and < h\ > are 
zero. The approximate estimate of the first order cor- 
rection to energy gives the nature of the phase near the 
boundary. However, since the estimate is only approxi- 
mate, at some regions of the parameter space, it is not 
reliable (as explained above). 
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APPENDIX A 
HIGHER DERIVATIVES FOR THE FIRST 
ORDER PERTURBATION CORRECTION 

We have derived earlier an expression 

E (i) = NT 2 sin O + 2irr - 3r0 o 
2-7T - O 

for the first order correction to the ground state energy 
of the Hamiltonian (see Eq. (3)) 



under perturbation by the Hamiltonian (see Eq. (4)) 



In this Appendix we shall show by the method of induc- 
tion, that for all values of 0o in the range < 0o < n, 
all higher derivatives of this quantity remain finite. 
For this purpose, we note that the first and second order 
derivatives are : 



aE« / 4tt 

= -AT — 7t -T- + 3 



and 



Or 



d 2 E^ 



(2tt - O ) 



-NT 



dr 2 \ (2tt — 0o ) 3 sin i/>o 

Here we have used the equality 
90o 2tt 



dr (2tt — 0o ) sin O 
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which can be obtained from Eq. (13). Thus, at least for 
m = 1 and 2 the derivatives can be expressed in the form 



d m E^ 



NY V k ^ 

h (27r-0o)*sinT^o 



(29) 



where m! ~ 3 m , and Qfj, 7,-, fcj are constants > 
with the restriction 7^ < a* (except for to = 1, where 
a, = 7j = 0). Differentiating Eq.(29) once, it can 
be easily verified that the form (29) should be valid 
for m + 1 also, and is hence valid for all values of m. 
Since the right hand side of Eq. (29) diverges only 
when sin 0o vanishes, we conclude that so long as (f>o is 
neither or tt the higher derivatives of E^ remain finite. 



APPENDIX B 
LONG RANGE CORRELATION FOR THE 
EXCITED STATES 



Appendix the expressions for the longitudinal two-spin 
correlation function 

C z (n) ^< s z jS * j+n > - < a* > 2 (33) 

(for n <C N) in the long-range limit n — > 00. 
Let us consider the quantity, 



a z (n) =< s z s z +n > 



(34) 



remembering that, 



lim o z (n) =< s z > 2 

Standard treatise [^tlH EH show that for = tt (no 
excited state), this quantity is of the following form. For 
A > 1, 



CT »=(1-- 



1/4 



Aexp(—an) 



Here we shall consider transverse Ising model with 
nearest neighbour interaction described by the Hamil- 
tonian 



A' 



w" = -]T[J^ +1 +r s *] 



(30) 



The exact solution 0, [H, [3| for this Hamiltonian tells 
us that the 2 N number of energy eigenvalues are 



(31) 



where A fc = ^/(A 2 + 2Acosfc + 1), A = J/T, ^ may be 0, 
±1 and k runs over N/2 equispaced values in the interval 
to vr. 

In the text we have come across (more than once) 
distribution of the following type : 



— 1 for k = to (f>. (unexcited) 
1 for k = 4> to 7r (excited) 



(32) 



This is an excited state and reduces to the ground state 
only for (f> = tt. This state is important for this work 
because it corresponds to the ground state after pertur- 
bation in the second perturbation scheme (Eq. (25)). It 
also corresponds to the first order correction in the first 
perturbation scheme (Eq. (11)). We shall present in this 



(A, a are constants depending on A.) This corresponds 
to ferromagnetic state. For A < 1, 

a z (n) = A' exp(— a'n) 

(A' , a' are constants depending on A.) This corresponds 
to paramagnetic state. There is a non-analyticity at the 
point A = 1 where the second derivative of ground-state 
energy diverges. 

Let us consider the correlation function C' z (n) when 
A ^ 1 and < <fi < tt (i.e. 0^0 and <f> ^ tt - some, 
but not all states are excited). We shall show elsewhere 
[l8| that this quantity can be calculated exactly using 
Szego's Theorem. The results obtained are as follows. 
The quantity linin^oo ct z (n) will be zero, so that C z (n) 
and cr z (n) are equal. Also, for A > 1, 



y/smq) \ A z 



1 \ cos(tt + <b)n 



(35) 



and for A < 1, 



C z (n) = 



0.590Vsirr; 



x ci +.v + 2\ <■<*<!>) ' (36) 



These expressions has been used in the text (Sec. IID, 
IIIB, IIIC) for the analysis of the phase diagram of ID 
transverse ANNNI model. 
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